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ABSTRACT 


The  classical  Nystrfcm  method  for  the  numerical  solution  of  Fredholm 
integral  equations  of  second  kind  consists  of  numerical  integration,  colloca- 
tion, and  interpolation.  The  approximate  solution  obtained  by  this  procedure 

is  shown  to  be  identical  to  the  solution  of  certain  finite-rank  integral 

/» 

equations  with  kernels  belonging  to  a specified  class  <K  and  thus  has 

n^ 

minimal  error  with  respect  to  approximation  of  the  original  equation  over 
this  class.  A computable  (but  in  general  nonoptimal)  error  bound  for  the 


Nystrom  approximate  solution  can  be  obtained  on  the  basis  of  how  well  a 

specific  finite-rank  integral  operator  with  kernel  in  (k  ^ approximate 

n 

integral  operator  in  the  Fredholm  equation  being  solved  numerically. 
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SIGNIFICANCE  AND  EXPLANATION 


Many  problems  in  applied  mathematics  (such  as  boundary-value  problems  for 
ordinary  and  partial  differential  equations)  can  be  solved  if  a function  x(s) 
satisfying  a Fredholm  integral  equation  of  second  kind, 

1 

(1)  x(s)  - A j K(s,t)x(t)dt  = y(s),  0 < s <_  1 , 

0 

can  be  found,  where  y(s)  and  K(s,t)  are  known.  In  most  cases,  x(s)  can- 
not be  found  explicitly,  so  numerical  methods  must  be  used.  E.  J.  Nystrom 
developed  a simple  method  of  this  type,  using  numerical  integration  to  replace 
the  integral  in  (1)  by  the  finite  sum 

n 

I K(s,t.)w  z , 
j=l  J J J 

where  z.  are  approximations  to  the  values  x(t.),  the  points  t, ,t  . ...,t 
3 3 1 2 n 

are  the  nodes  and  w, ,w  . ...,w  the  weights  of  the  rule  of  numerical  integra- 

12  n 

tion  (Simpson,  Gauss,  etc.)  being  used.  In  order  to  determine  the  approximate 

values  z, , z„,...,z  , collocation  at  the  points  t, ,t„,...,t  is  used  to 
1 2 n 1 2 n 

obtain  the  system  of  equations 

n 

(2)  z - A l K(t.,t.)w  z.  = y(t.),  i = 1,2, ...,n  , 

j=l  1 J J J 

to  be  solved  for  these  values.  Finally,  interpolation  by  the  formula 

n 

(3)  z(s)  = y (s)  + X l K(s,t.)w.z. 

j=l  3^3 

gives  an  approximate  solution  of  equation  (1)  on  the  entire  interval  0 <_  s <_  1. 
This  method  is  simple  in  concept  and  easy  to  implement  on  electronic  computers. 

The  question  remains,  how  accurate  are  the  results?  In  fact,  the  Nystrom  method 
is  known  to  produce  good  approximate  solutions  in  many  applications,  but  why  is 
this  so?  In  this  report,  a number  of  ways  of  analyzing  the  accuracy  of  the 
Nystrom  method  are  compared,  and  it  is  shown  that  error  can  be  estimated  by  a 
simple  computational  procedure  based  on  approximation  of  K(s,t)  by  kernels 
of  the  form 

n 

(4)  K ( s , t ) - l K(s,t.)w  v (t)  . 

n j*»l  J J J 

It  is  also  shown  that  the  error  of  the  Nystrom  method  depends  on  the  best  possible 
approximation  of  K(s,t)  by  kernels  Kn(s,t)  from  a given  class,  and  thus  the 
observed  accuracy  of  this  method  has  a theoretical  explanation. 

The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


For  if  K is  a linear  integral  operator  in  the  apece  CIO, 11  with  operator  nor*  of 


1 


(7.3) 


oiiginal  task  of  solving  the  linear  mteoral  equation  (1.1). 
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